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The unambiguous detection of Galactic dark matter annihilation would un- 
ravel one of the most outstanding puzzles in particle physics and cosmology. 
Recent observations have motivated models in which the annihilation rate is 
boosted by the Sommerfeld effect, a non-perturbative enhancement arising 
from a long range attractive force. Here we apply the Sommerfeld correc- 
tion to Via Lactea II, a high resolution TV-body simulation of a Milky- Way-size 
galaxy, to investigate the phase-space structure of the Galactic halo. We show 
that the annihilation luminosity from kinematically cold substructure can be 
enhanced by orders of magnitude relative to previous calculations, leading to 
the prediction of 7-ray fluxes from up to hundreds of dark clumps that should 
be detectable by the Fermi satellite. 

In the standard cold dark matter (CDM) paradigm of structure formation, a weakly inter- 
acting massive particle (WIMP) of mass m x ~ 100 GeV-10 TeV ceases to annihilate when the 
universe cools to a temperature of Tf ~ m x /20, about one nano-second after the Big Bang. A 
thermally-averaged cross-section at freeze-out of (av)o ~ 3 x 10~ 26 cm 3 s _1 results in a relic 
abundance consistent with observations (/). Perturbations in the dark matter density are ampli- 
fied by gravity after the universe becomes matter dominated, around ten thousand years after 
the Big Bang: the smallest structures ("halos") collapse earlier when the universe is very dense 
and merge to form larger and larger systems over time. Today, galaxies like our own Milky 
Way are embedded in massive, extended halos of dark matter that are very lumpy, teeming with 
self-bound substructure ("subhalos") that survived this hierarchically assembly process (2-4). 
Indirect detection of high energy antiparticles and 7-rays from dark matter halos provides a po- 
tential "smoking gun" signature of WIMP annihilation (5). The usual assumption that WIMP 
annihilation proceeds at a rate that does not depend, in the non-relativistic «/c C 1 limit, on 
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the particle relative velocities implies that the primary astrophysical quantity determining the 
annihilation luminosity today is the local density squared. WIMP annihilations still occur in the 
cores of individual substructures, but with fluxes that are expected to be dauntingly small. The 
latest calculations show that only a handful of the most massive Galactic subhalos may, in the 
best case, be detectable in 7-rays by the Fermi satellite (6, 7). 

The Sommerfeld enhancement, a velocity-dependent mechanism that boosts the dark matter 
annihilation cross-section over the standard (av) value (8-11), may provide an explanation for 
the experimental results of the PAMELA satellite reporting an increasing positron fraction in 
the local cosmic ray flux at energies between 10 and 100 GeV (12), as well as for the surpris- 
ingly large total electron and positron flux measured by the ATIC and PPB-BETS balloon-borne 
experiments (13, 14). Very recent Fermi (15) and H.E.S.S. (16) data appear to be inconsistent 
with the ATIC and PPB-BETS measurements, but still exhibit departures with respect to stan- 
dard expectations from cosmic ray propagation models. Although conventional astrophysical 
sources of high energy cosmic rays, such as nearby pulsars or supernova remnants, may provide 
a viable explanation (17-19), the possibility of Galactic DM annihilation as a source remains in- 
triguing (20-22). In this case, cross-sections a few orders of magnitude above what is expected 
for a thermal WIMP are required (23). 

The Sommerfeld non-perturbative increase in the annihilation cross-section at low velocities 
is the result of a generic attractive force between the incident dark matter particles that effec- 
tively focuses incident plane- wave wavefunctions. The force carrier may be the W or Z boson 
of the weak interaction (10), = 80 — 90 GeV/c 2 , or a lighter boson, ~ GeV/c 2 , mediating 
a new interaction in the dark sector (11, 24). Upon introduction of a force with coupling strength 
a, the annihilation cross-section is shifted to (av) = S(av) , where the Sommerfeld correction 
S disappears (S = 1) in the limit v/c — > 1 (thus leaving unchanged the weak scale annihilation 
cross-section during WIMP freeze-out in the early universe). When v/c <C a, S « irac/v 
enhancement), but it levels off to S max « 6am x /m^ at v /c « 0.5m ( j ) /m x because ofthe 
finite range of the interaction. For specific parameter combinations, i.e. when m x /m ( f > m n 2 / a 
where n is an integer, the (Yukawa) potential develops bound states, and these give rise to large, 
resonant cross-section enhancements where S grows approximately as 1/v 2 before saturating 
(see Supporting Online Material). 

The Sommerfeld effect connects dynamically the dark and the astrophysics sectors. Be- 
cause the typical velocities of dark matter particles in the Milky Way today are of the order of 
v/c ~ 10~ 3 , the resulting boost in the annihilation rate may provide an explanation to the puz- 
zling Galactic signals. Compared to particles in the smooth halo component, the Sommerfeld 
correction preferentially enhances the annihilation luminosity of cold, lower velocity dispersion 
substructure, as emphasized previously by (10, 25, 26). Detailed knowledge of the full phase- 
space density of dark matter particles in the Milky Way is thus necessary to reliably compute 
the expected signals. Here we use the Via Lactea II cosmological simulation, a high preci- 
sion calculation of the assembly of the Galactic CDM halo, for a systematic investigation of 
the impact of Sommerfeld-corrected models on present and future indirect dark matter detec- 
tion efforts. Via Lactea II employs just over one billion 4, 100 M Q particles to follow, with a 
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Figure 1: A: The distribution of velocity dispersion for Via Lactea //particles within 400 kpc. The dispersions 
are calculated from the nearest 32 neighbors of each particle. B: Sommerfeld enhancement factor as a function 
of velocity for four representative models, exhibiting S ~ 1/v (cyan and orange curves) and ~ 1/v 2 behavior 
(magenta and brown), and high (cyan and magenta) versus low (orange and brown) saturation velocities. C and D: 
The corresponding distributions of S-factors for Via Lactea //particles. 

force resolution of 40 pc, the formation of a 1.9 x 10 12 M Milky-Way size halo and its sub- 
structure from redshift z = 104 to the present (27-29). (Fig. 1 A) The smooth halo particles, 
whose velocity dispersions are set by the global potential, typically have three-dimensional ve- 
locity dispersion o > 100 kins" 1 . Particles in self-bound subhalos dominate at lower velocity 
dispersions. The total mass fraction of particles with a < 5 kms" 1 is 1%. We calculated 
the Sommerfeld enhancement factors S on a particle-by-particle basis by averaging S(v) over a 
Maxwell-Boltzmann distribution of relative velocities with one-dimensional velocity dispersion 
given by y / 2/3 a (see the Supporting Online Materials for details) (Fig. 1 C, D). 

The large Sommerfeld boost expected for v/c ~ 10~ 4 — 10~ 5 make cold subhalos more 
promising sources of annihilation 7-rays than the higher density but much hotter region around 
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Figure 2: All-sky maps (in a Mollweide projection) of the Sommerfeld-enhanced annihilation surface brightness 
(Jios from all Via Lactea II dark matter particles within 400 kpc. The observer is located at 8 kpc from 

the halo center along the host halo's intermediate principal axis. A: No Sommerfeld enhancement. B: S ~ 1/v, 
saturated at ~ 1 km s — 1 . C: S ~ 1/v 2 saturated at ~ 5 km s _1 . The maps have been normalized to give the same 
total smooth host halo flux. 

the Galactic Center (Fig. 2). In Sommerfeld-enhanced models, substructures are much more 
clearly visible, and can even outshine the Galactic Center when the cross-section is close to 
resonance and saturates at low velocities. Furthermore, baryonic processes will tend to heat up 
the Galactic Center and dim its Sommerfeld boost, and thereby increase the relative detectability 
of subhalos. Dark matter halos are not isothermal and have smaller velocity dispersions in the 
center (see Supporting Online Material). In addition to an overall increase in the annihilation 
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Figure 3: Annihilation rate maps (projections of p 2 S out to the tidal radius) of one of the most massive (M ~ 

2 x 10 9 M Q ) subhalos in Via Lactea II, for the same models as in Fig. 2. The images have not been normalized. 
Compared to the S = 1 case (A), the total luminosity is 2,200 and 160,000 higher for S ~ 1 /v (B) and S ~ l/v 2 
(C), respectively. 

rate, this "temperature inversion" leads to a relative brightening of the center at the expense of 
the diffuse flux from the surrounding region (Fig. 3). The subhalo exhibits its own population 
of subclumps, also Sommerfeld-enhanced. 

To address quantitatively the detectability of Sommerfeld-enhanced subhalos by the Fermi 
Space Telescope, we have converted the annihilation flux calculated from our simulation (30) 
into a predicted 7-ray flux and compared it to the expected backgrounds. We investigated two 
different classes of particle physics models (Table 1): 

i) those motivated by (70), in which the force carrier is the conventional weak force gauge 
boson, the W or Z particle, and the mass of the dark matter particle is > 4 TeV. We have 
chosen four representative values of m x and a, which lie increasingly close to a S ~ l/v 2 
resonance. In these models the main source of 7-rays is the decay of neutral pions that 
are produced in the hadronization of the annihilation products; 

ii) and those in which the annihilation is mediated by a new dark sector force carrier (f) (11). 
The choice of parameters (m x , m^, a) follows Meade, Papucci, & Volansky (MPV) (31) 
and satisfies constraints from recent H.E.S.S. measurements of the Galactic Center (32) 
and the Galactic Ridge (33), as well as the PAMELA measurement of the local positron 
fraction above 10 GeV (12) and the ATIC (13) and PPB-BETS (14) measurements of the 
total (e + + e~) flux above 100 GeV. Models MPV-1 incorporate all three constraints and 
models MPV-2 only the H.E.S.S. and PAMELA data. We considered models away from 
(a) and close to (b) resonance, thereby covering both S ~ l/v and ~ l/v 2 behavior. The 
data favor a light force carrier, ps 200 MeV, and the 7-rays originate then as final state 
radiation (internal bremsstrahlung) accompanying the decay of the 0's into e + e~ pairs. 

The magnitude of the relativistic cross section was fixed to the standard value of (crv)o = 

3 x 10~ 26 cm 3 s _1 . 7-ray spectra are shown in Figure S2 in the Supporting Online Material. 
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Table 1 : Summary of the models used to assess subhalo detectability with Fermi. Particle physics parameters 
are: m x , the mass of the dark matter particle, m^, the mass of the force carrier, and a, the coupling constant. In 
the two right-most columns we give S max , the maximum Sommerfeld enhancement obtained, and the saturation 
velocity v sa ±, defined as the velocity at which S reaches 90% of S m&x . 



Model 


m x 


m 


a x 100 


<? 

'-'max 


^sat 




(TeV) 


(GeV) 






(km s" 1 ) 


LS-1 


4.30 


90 


3.307 


1,500 


80 


LS-2 


4.45 


90 


3.297 


12,000 


28 


LS-3 


4.50 


90 


3.288 


70,000 


12 


LS-4 


4.55 


90 


3.281 


430,000 


4.7 


MPV-la 


1.0 


0.2 


4.000 


3,000 


7.4 


MPV-lb 


1.0 


0.2 


3.739 


16,000 


2.4 


MPV-2a 


0.25 


0.2 


4.000 


480 


40 


MPV-2b 


0.25 


0.2 


4.500 


40,000 


3.3 



We determined the Fermi detection significance by summing the annihilation photons from 
all the pixels in our all-sky maps covering a given subhalo, and compared this to the square root 
of the number of background photons from the same area. We counted a subhalo as "detectable" 
if it had a total signal-to-noise greater than 5 (Table [2]). The numbers are quite large, implying 
that individual subhalos should easily be detected by Fermi if Sommerfeld enhancements are 
important. Even in the most conservative cases (MPV-la and MPV-2a) around ten or more sub- 
halos should be discovered after 5 years of observation. In fact, on the basis of all models con- 
sidered here it is predicted that Fermi should be able to accumulate enough flux in its first year 
of observations to detect several dark matter subhalos at more than 5a significance, a prediction 
that will soon be tested, and that may open up the door to studies of non-gravitational dark 
matter interactions and new particle physics. The central brightening discussed above results 
in a smaller angular extent of a given subhalo's detectable region: the stronger the Sommer- 
feld enhancement the fewer pixels exceed the detection threshold. Nevertheless, for all models 
considered here the majority of detectable subhalos would be resolved sources for Fermi. 

Another question of interest is whether Sommerfeld-corrected substructure would lead to 
a significant boost in the local production of high energy positrons, arising from dark matter 
annihilation in subhalos within a diffusion region of a few thousand parsecs from Earth, as 
well as of antiprotons within a correspondingly larger diffusion region. The local dark matter 
distribution at the Sun's location appears quite smooth in the highest-resolution numerical sim- 
ulations to date (27, 29, 34). Tidal forces efficiently strip matter from subhalos passing close 
to the Galactic Center and often completely destroy them. Further substructure depletion may 
be expected from interactions with the stellar disk and bulge. In the Via Lactea II simulation, 
the mean number of > 10 5 M Q subhalos within 1 kpc of the Sun is only 0.04, and one must 
reach three times farther to find one clump on average. Without the Sommerfeld effect, this 
dearth of nearby substructure leads to a local annihilation boost of less than 1%, and at most 
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Table 2: Detectable Subhalos. The number of subhalos that would be detected with > 5a significance by Fermi 
after 1, 2, 5, and 10 years in orbit, for different Sommerfeld-enhanced dark matter particle models. In the two 
right-most columns we give the median distance and mass of the detectable clumps after 5 years in orbit. 



Model 


1 yr 


2 yr 


5 yr 


10 yr 


D 


Msub 












(kpc) 


(M ) 


LS-1 


12 


19 


29 


38 


24 


1.4 x 10 7 


LS-2 


72 


99 


167 


244 


42 


9.5 x 10 6 


LS-3 


225 


311 


457 


583 


56 


6.2 x 10 6 


LS-4 


410 


528 


730 


919 


66 


4.9 x 10 6 


MPV-la 


5 


7 


12 


15 


16 


9.8 x 10 7 


MPV-lb 


9 


14 


25 


36 


25 


4.4 x 10 6 


MPV-2a 


12 


18 


29 


38 


24 


1.4 x 10 7 


MPV-2b 


187 


254 


397 


518 


55 


4.5 x 10 6 



20% in the rare case of a nearby clump, as found in a statistical approach (35). The picture 
changes with Sommerfeld enhancement. The low velocity dispersion of cold substructure leads 
to a greatly increased luminosity compared to the hotter smooth component. For typical 1/v 
models, subhalos resolved in our simulation within 2 kpc contribute on average about half as 
much luminosity as the smooth component, and up to 5 times as much in rare cases. If the 
Sommerfeld enhancement is resonant (S ~ 1/v 2 ), then these subhalos dominate by a factor of 
20 on average and by as much as 200 in rare cases. 
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Figure S 1 : A: The Sommerfeld enhancement factor S as a function of m-^/m^ (the mass ratio of the dark matter 
particle to the force carrier particle) at a fixed coupling (a = 0.30) for different velocities. B and C: S as a function 
of velocity for the models for which we calculate subhalo detectability with Fermi. The parameters of the models 
are given in Table 1 in the main text. 



A Sommerfeld enhancement to the annihilation cross-section arises when the dark matter 
(DM) particle is heavy compared to the gauge boson mediating the interaction: m x > m^/a, 
where a = A 2 /47r, and A is the coupling between the dark matter particle \ an d the force 
carrier (p. The magnitude of the enhancement can be determined by solving the two-body radial 
Schrodinger equation, as shown in (1, 2). Defining the additional dimensionless parameter 



P* 



1 00710 



(1) 



one can distinguish three regimes for the dependence of the enhancement S as a function of the 
velocity (2): 

i) for large velocities, [3 = v/c ^> a, there is no Sommerfeld enhancement, S ~ 1. This 
ensures that the relic abundance of dark matter is not affected, since the velocities were 
close to relativistic at freeze-out; 

ii) at intermediate velocities, <C <C a, the Sommerfeld enhancement follows a 1/v 
behavior, S w not/ (3; 
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iii) at low velocities, (3 <C (3*, resonances appear for certain values of m x , due to the pres- 
ence of bound states. In this case S grows as 1/v 2 before saturating at a velocity that 
depends on how close to resonance m x lies. In the non-resonant case the Sommerfeld en- 
hancement saturates when the deBroglie wavelength of the particle ~ (m x v)~ l becomes 
comparable to the range of the interaction ~ mT 1 ), i.e. when v ~ m^jm x . 



Figure [ST] shows the rich behavior of the Sommerfeld enhancement, obtained by numerical 
integration of the Schrodinger equation with a Yukawa potential. 

The magnitude of the Sommerfeld enhancement at a given location depends on the details 
of the distribution of relative velocities of the DM particles. Accurately determining the phase- 
space structure as a function of position is a notoriously difficult task even with the highest 
resolution N-body simulations. Some recent investigations based on cosmological N-body sim- 
ulations have found significant structure in the coarse-grained velocity distribution of the host 
halo and departures from a simple, single- "temperature" Maxwell-Boltzmann distribution (3,4). 
These features are remnants of the accretion history of the host halo due to incomplete phase 
mixing. The centers of subhalos, however, are likely more well-mixed and may show less sig- 
nificant departures from a Maxwellian. At any rate, a full characterization of the phase-space 
distribution at all locations in our simulation is beyond the scope of this paper, and hence we 
make the simplifying assumption that the distribution of relative velocities is of the Maxwell- 
Boltzmann form 



f(Vrel] = 4:71 



(2trx Mi1D ) 3 / 2 



v rel exp 



1 / ^rel 

2 \Cu,lD. 



(2) 



with a one-dimensional velocity dispersion 

kT 



a 



kT 



m 



2 (X, 



m,lD) 



(3) 



where \i = mim 2 /(mi + m 2 ) = m/2 is the reduced mass of a two particle system. From 
the simulated particles we determine a three-dimensional velocity dispersion, a 2 = a, 



m,3D 



3 ^m, ID 



3/2o"^ 1D (for systems with zero velocity anisotropy). The Maxwell-Boltzmann- 
weighted Sommerfeld enhancement is then given by 



S(a) 



f{v-^j2/ia) S(v) dv. 



(4) 



S2 Gamma-rays from dark matter annihilation 

Since the dark matter particle is neutral it does not couple directly to the electromagnetic field, 
and hence annihilations straight into two monochromatic photons (or a photon and a Z boson) 
are typically strongly suppressed. Nevertheless 7-rays can be a significant by-product of dark 
matter annihilations, since they can arise either from the decay of neutral pions produced in 
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Figure S2: The 7-ray spectrum per annihilation for the models under consideration here. In the LS models 
(XX ~ * ZZ or W + W~) the 7-rays come from the decay of pions produced in the decay of the bosons. The dotted 
line indicates the internal bremsstrahlung contribution in the case of W + W~ (5). In the MPV models (xX ~* ( t )( t ) ) 
the 7-rays originate as final state radiation associated with the decay of the tj> carriers into e + e~ pairs. 

the hadronization of the annihilation products, or through internal bremsstrahlung associated 
with annihilations into charged particles, or from interactions of energetic leptons with the 
surrounding interstellar photons (inverse Compton scattering). We do not consider the latter 
process here, since we are focusing on the annihilation signal from dark matter subhalos which 
are unlikely to harbor a sufficiently high stellar radiation field. 

The 7-ray spectra per annihilation that we use in our detectability calculation are shown in 
Figure [S2] In the Lattanzi & Silk models the annihilation results in two neutral Z bosons or 
a pair of W + and W~ bosons, and the dominant source of 7-rays is neutral pion decay. For 
m x = 4.5 TeV, every annihilation results in ~ 26 photons with energies between 3 and 300 
GeV. In the MPV models, the mass of the <p particle is so low (by design), that only decays 
into e + e~ pairs are kinematically allowed, and we must rely on final state radiation (internal 
bremsstrahlung) for the 7-ray signal. This results in fewer 3-300 GeV 7-rays per annihilation 
(0.39 for m x = 1 TeV, 0.30 for m x = 250 GeV), but this is partially compensated by the smaller 
dark matter particle mass and hence higher number density at fixed mass density. 

S3 The velocity structure of the Via Lactea II host and its 
subhalos 

In Figure [S3] we present radial profiles of the density p, velocity dispersion a, and velocity 
anisotropy parameter j3 for the smooth host halo and averaged over the 100 most massive sub- 
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Figure S3: Left panel: The radial dependence of the density (A), velocity dispersion (B), and the velocity 

anisotropy (3 = 1— ^ e $ * (C) for the smooth host halo component. The solid red line shows the values 
calculated from all particles in spherical shells. For the p and a profiles we also show the median (dark green 
dashed line) and the 68% region (light green shaded region) of the particle -by-particle quantities determined from 
the nearest 32 neighbors. The vertical dashed line indicates our estimate for the convergence radius of the density 
profile (380 pc). The shape of the a profile implies that Sommerfeld enhancement will preferentially brighten the 
very central region of the host halo at the expense of the surrounding region. The anisotropy profile clearly shows 
that the host halo is not isotropic, with a slight radial anisotropy persisting down to the convergence radius. Right 
panel (D-F): The same quantities averaged over the 100 most massive subhalos. The bullets indicate the median 
and the error bars the 68% scatter around the median. In D and E we also plot the best-fitting NFW (solid) and 
Einasto (a fixed at 0. 17, dotted) profiles. Note that the subhalo profiles should not be considered converged below 

~ 0.1 r V max- 



halos. We determine these profiles by first binning all particles into equally spaced logarithmic 
radial shells, and then calculating 



P(r) = (5) 



Anr 2 dr 

-2/„\ _ //„. /„. \ \2 



°#r) = ((vj - (vjhYh, (6) 
a\r) = a 2 x (r) + cr 2 y (r) + a 2 z (r), (7) 



1 crf(r) + al(r) 

2 erf r 



where the sum and averages (denoted by are over all particles in a given spherical shell. 
These profiles are indicated by the solid red line for the host halo in the left panels of Figure IS3l 
For the p(r) and a(r) profiles we also show the median and 68% interval of the distribution 
of particle density pi and velocity dispersion <7i, both calculated from the 32 nearest neigh- 
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boring particles. The two estimates agree quite well with each other, with the slightly higher 
(lower) median p { (c^) indicating a negative (positive) skew of the distribution, presumably due 
to spherical averaging of a triaxial mass distribution. Since 32 particles are not sufficient for a 
good estimate of (3, we don't show the distributions for the velocity anisotropy. 

The density profile of the Via Lactea II host has been discussed in (6) and we present it 
here merely for completeness. Down to our convergence radius of 380 pc it is well fit by a 
generalized NFW profile, 

p s 2 3 ~ 7 

= ( I Vy/i i / (9) 

(r/r s p(l + r/r s ) 6 7 



with a central slope of 7 = 1.2, but an Einasto profile 

p(r) = p s exp 



-f ((r/r s y 



(10) 



with a = 0.167 fits almost as well. The velocity dispersion profile exhibits the central "tem- 
perature inversion" typical of cold dark matter halos: er(r) peaks at about the scale radius and 
decreases with decreasing radius (7-9). The /3(r) profile shows the well established trend of a 
considerable amount of radial anisotropy in the outskirts of the halo and decreasing towards the 
center (7-9). We find that even at the convergence radius a slight amount of radial anisotropy 
remains. 

In the right panels of Figure[S3] we show the p, a, and (3 profiles averaged over the 100 most 
massive subhalos in our simulations. We scaled the radius of each subhalo by its rv max and 
calculated radial profiles using 30 equally spaced logarithmic bins from r/r Vm ax = 0.01 to 2. 
In each radial bin we then determined the median and 68% region of the distribution of values 
over all 100 subhalos, rejecting bins containing less than 100 particles. We only plotted bins 
containing values from more than 10 subhalos. It is difficult to estimate a convergence radius 
for these subhalos; the host halo convergence radius of 380 pc corresponds to (0.05 - 0.5) r Vm ax 
for these 100 subhalos, so these average profiles should not be considered converged below 

~ O.lrvmax- 

The median density profile nicely follows the anticipated NFW-like profile. We have over- 
plotted the best-fitting NFW (solid line) and Einasto (with fixed a = 0.17, dotted line) profiles. 
Clearly the resolution is not good enough to allow a quantitative analysis of the asymptotic cen- 
tral slope of the density profile, but it remains cuspy as far down as we can resolve. The velocity 
dispersion profile looks qualitatively the same as the host halo's, with a peak around 0.25r Vm ax 
and decreasing towards smaller radii. Not surprisingly, subhalos are not isothermal. The strong 
increase in the 84 th percentile of a at large radii is due to contamination by host halo particles 
which artificially inflate a. Although the (3 profile is quite noisy, it is clear that subhalos too 
exhibit a slight radial anisotropy even down at the smallest radii that we can access. 

An important caveat to these findings is that our simulations completely neglect the effects 
of baryons on the DM distribution. Gas cooling, star formation, supernova feedback, and stellar 
dynamical processes might alter both the DM density and velocity dispersion profiles. In fact, 
not even the sign of these effects is clear at the moment: adiabatic contraction generally leads 
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to a steepening of the central density profile (10, 11), but dynamical friction acting on baryonic 
condensations tends to remove the central cusp (12, 13). The velocity dispersion of DM par- 
ticles, on the other hand, is more likely to increase in regions affected by baryonic processes. 
These complications will be most important for the Galactic Center. The high mass-to-light ra- 
tios observed in Galactic dwarf satellites (14) indicate that they are completely DM dominated 
and hence likely much less affected by baryonic physics. Interactions with the Milky Way's 
stellar and gaseous disk probably strip a significant fraction of mass from some DM subhalos, 
but the dense central regions responsible for most of the annihilation luminosity are relatively 
well protected. 



S4 Detectability Calculation 

We calculated the annihilation flux directly from our simulations following the procedure de- 
tailed in (16), with one important modification to account for the Sommerfeld enhancement. 
The intensity in a given direction (6, <p) is now given by 

l(9,cj))=Q p 2 S(a;m x ,m (j} ,a)d£, (11) 

Jlos 

where Q contains most of the particle physics dependence (the particle mass, the high velocity 
annihilation cross section, and the 7-ray spectrum per annihilation event) and S(a; m x , m^, a) 
is the Maxwell-Boltzmann- weighted Sommerfeld enhancement factor at a velocity dispersion a 
for a given particle physics model. For discrete particles, with masses mi, distances dj, densities 
Pi, and velocity dispersions <7;, this integral becomes a discrete sum over all particles in a given 
map pixel, 

\- piS(cTi)mi 
i 4< • (12) 

We determined pi and <7j from the nearest 32 neighbors of the ith particle, but have checked 
that our results do not change significantly for 64 neighbors. We have implemented several 
additional improvements over (16): (a) we corrected a calculation error in the conversion from 
simulation fluxes to gamma-ray counts (the subhalo fluxes were a factor of two too large in 
(16)), (b) we use a 15% lower effective area (~ 7,300 cm 2 ), as suggested by the performance 
of the LAT instrument measured in orbit (17), and (c) we switched to the HEALPixQ equal 
area pixelization scheme, setting N side = 512, which corresponds to a solid angle per pixel of 
Afl = 4 x 10~ 6 sterad, comparable to the angular resolution of Fermfs LAT detector above 3 
GeV. We assumed a LAT exposure time equal to 0. 153 of the time in orbit, a combination of the 
~ 47r/5 sr field of view, 90% trigger live time, and 15% data acquisition outage during South 
Atlantic Anomaly passages (17). Our analysis was restricted to one fiducial observer located 
at 8 kpc from the host halo center along the intermediate axis of its density ellipsoid, and we 

'http : / /healpix. jpl .nasa . gov/ 
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refer the reader to (16) for a discussion of the signal variance arising from different observer 
locations. 

Due to the finite resolution of our simulation the centers of all our halos are artificially heated 
and less dense than they would be at higher resolution. This results in central brightnesses that 
are lower than would be expected for an NFW or Einasto density profile. We have corrected the 
central flux from the host halo and all subhalos using an analytical extrapolation of the density 
and velocity dispersion profiles. We considered both an NFW (7 = 1) profile and an Einasto 
profile with a = 0.17, which has been shown to fit Galactic-scale dark matter halos (18). These 
analytical profiles are matched to the measured values of \/ max and r Vm ax, which are robustly 
determined for subhalos with more than 200 particles (M > 8 x 10 5 M ). The relations between 
(Knax, ?vmax) and (p s , r s ) are 

r V max = fr f s V^ax = f V ^Gp^, (13) 

with f r ~ 2.163 (2.212) and f v « 0.865 (0.897) for the NFW (Einasto) profile. We use the 
spherical Jeans equation to solve for the corresponding velocity dispersion profile, assuming 
(3 = 0. For a halo at a given distance we can then solve for a Sommerf eld-enhanced surface 
brightness profile as a function of angle from the halo center, average it over the angular reso- 
lution of our maps (Af2 = 4 x 10~ 6 sterad), and use this to correct our simulated maps. For 
the host halo we only correct pixels within ~ 2.7° from the center, corresponding to the den- 
sity profile convergence radius of 380 pc. For the subhalos we ensure that all pixels within the 
projected scale radius r s have a surface brightness at least as high as the expectation from the 
analytical extrapolation. These correction factors are typically not very large: over all subhalos 
and all Sommerfeld models, the median (root mean square) correction factor for the central 
pixel is 2.2 (5.0). 

Note that we neglect the possible enhancement of a subhalo's luminosity arising from ad- 
ditional clumpy substructure below our simulation's resolution limit. The magnitude of this 
so-called "substructure boost factor" depends on uncertain extrapolations of the abundance, 
distribution, and internal properties of the low mass subhalos, and will not be significantly 
increased by the Sommerfeld effect due to its saturation at low velocities. While this boost 
typically doesn't affect the central surface brightness very much, it may somewhat increase the 
angular extent of a given subhalo's signal. 

The annihilation signal from individual subhalos must compete with a number of diffuse 
7-ray backgrounds, of both astrophysical and DM annihilation origin. At low Galactic lati- 
tudes the dominant astrophysical background arises arises from the interaction of high energy 
cosmic rays with interstellar gas (pion decay and bremsstrahlung) and radiation fields (inverse 
Compton). This background has been measured by the EGRET instrument aboard the Comp- 
ton satellite at 0.5 degree resolution out to ~ 30 GeV (19). An improved measurement of the 
spectral and angular properties of this background is one of the goals of the Fermi mission, and 
preliminary data at intermediate Galactic latitudes have already been presented by the Fermi 
collaboration (20). Here we employ a theoretical model of this background (the GALPROP 
"conventional" model (21, 22)), which matches the EGRET and preliminary Fermi measure- 
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Figure S4: Cumulative distribution of the angular size of the detectable subhalos. We plot the fraction of 
S/N> 5 subhalos with more than N p i x pixels exceeding the Fermi detection threshold after 5 years in orbit, for the 
LS models in panel A, and for the MPV models in panel B. 

merits with sufficient accuracy for our purposes. At high Galactic latitudes (|6| <: 20°), an 
isotropic extragalactic 7-ray background from unresolved blazars may become dominant. We 
include such a background with an intensity and spectrum as measured by EGRET (23). Note 
that this is probably an overestimate, since Fermi will likely resolve some fraction of this back- 
ground into point sources. We do not include the contributions from astrophysical foreground 
and background point sources, but we have checked that none of the Sommerfeld-enhanced sub- 
halos would have been detected by EGRET, assuming a point-source sensitivity of 2 x 10~ 7 7 
cm 2 s- 1 from 0.1 to 30 GeV (24). 

In addition to these astrophysical backgrounds, individually detectable subhalos must com- 
pete with the diffuse background from DM annihilation in the smooth host halo and from the 
population of individually undetectable subhalos (see Section |S5l) (16, 25). In models without 
Sommerfeld enhancement and for S ~ 1/v models, these DM annihilation backgrounds are 
negligible compared to the astrophysical backgrounds (although they themselves constitute a 
signal worth searching for). In resonant S ~ 1/v 2 models with low saturation velocity (e.g. 
LS-3, LS-4, MPV- lb, MPV- 2b), however, the unresolved subhalo flux becomes the dominant 
background and must be accounted for. For completeness we have included all four components 
in the detectability calculation in all cases. 

Subhalo detectability for Fermi is assessed as follows. First we calculate a signal-to-noise 
ratio (S/N) per pixel by dividing the number of source photons arising in the subhalo map by 
the square root of the number of photons in the background map. Next we select all pixels 
with S/N > 1 and identify contiguous regions, which we associate with individual subhalos 
based on proximity of the brightest pixel with a subhalo center. If more than one subhalo center 
coincides with a given brightest pixel, we pick the subhalo with the larger expected surface 
brightness (L/r 2 ~ KiLx/ r vmax)- F° r eacn °f the se contiguous regions we then calculate a 



17 



subhalo detection significance 

N 

S/N = 4=^, (14) 

where N s and N b are the total number of source and background 7-rays over the contiguous 
region. This definition is a good proxy for detection significance under the assumption that an 
estimate of the background can be subtracted out and only Poisson fluctuations remain. Note 
that it is not the uncertainty of the flux itself (which would be N s / y/N s + Nb), but instead an 
estimate of the significance of having detected a departure from a smooth background. 

In Figure [S4] we show the cumulative distribution of the angular size (the number of pixels 
exceeding the Fermi detection threshold after 5 years in orbit) of the detectable subhalos. Al- 
though models with stronger Sommerfeld enhancement result in smaller detectable regions, the 
majority of all subhalos would still be resolved sources for Fermi. This effect raises the possi- 
bility of future observations being able to discriminate between different Sommerfeld-enhanced 
models. 



S5 Diffuse Flux from Unresolved Subhalos 

For a typical CDM power spectrum of density fluctuations one would expect dark matter clumps 
on scales all the way down to a cutoff set by collisional damping and free streaming in the 
early universe (26, 27). For WIMP dark matter, typical cut-off masses are m = 1CT 12 to 
10" 4 M Q (28, 29), some 10 to 20 orders of magnitude below Via Lactea IPs mass resolution. 
In this case the Galactic dark matter halo might host an enormous number (10 10 — 10 22 , see 
(16)) of small mass subhalos, whose combined annihilation signal could result in a sizeable 
7-ray background. If the Sommerfeld enhancement didn't saturate at a finite velocity, this 
background would easily outshine any other Galactic 7-ray signal. Even with saturation one 
must ask whether the Sommerfeld-enhanced annihilation background from this population of 
unresolved subhalos would outshine individual subhalos. 

In order to address this question, we have extended the analytical model of (16) to allow 
for Sommerfeld enhancement. The overall intensity of this background depends sensitively on 
a number of uncertain parameters governing the subhalo population as a whole (slope and nor- 
malization of the mass function, their dependence on distance to the host center) and the mass 
dependent properties of individual subhalos (density and velocity dispersion profiles, concentra- 
tions). Although the model is calibrated to numerical simulations at the high mass end, it relies 
on an extrapolation over many orders of magnitude in mass below the simulation's resolution 
limit that is very uncertain. 

Our model employs a radially anti-biased subhalo mass function 

**£ri = 1.05 x 10-M5' kpc- (4Vr (l + — ^V' <'5) 
dM P V 10 M ©/ V 18.5 kpcy V 
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Figure S5: The annihilation intensity as a function of angle ip from the Galactic Center, for the host halo (red), 
individual subhalos (blue), and unresolved subhalos (magenta dashed). A: No Sommerfeld enhancement. B: 
Model MPV-2a (S ~ 1/v, high u sat ). C: Model LS-4 (S ~ l/v 2 , low u sat ). For the host halo and unresolved 
components we plot the mean, for the individual subhalo profile the maximum intensity over all pixels in a given 
ip bin. Individual subhalos outshine the diffuse unresolved subhalo background, even in the strongly Sommerfeld- 
enhanced case. 



with a low mass cutoff of m = 10 6 M Q , and assumes an Einasto density and velocity dis- 
persion profile with a concentration-mass relation according to (30) and a radial dependence 
of 



r 



-0.286 



^>-»(M) J^l (16) 



We refer the reader to the Appendix of (16) for details about the calculation. In Figure 
we show the resulting azimuthally averaged intensity as a function of angle from the Galactic 
Center, and compare it to the smooth host halo signal and the peaks due to individual subhalos. 
While the unresolved subhalo background is brighter than the smooth halo component every- 
where but in the very center, individual subhalos have higher central surface brightnesses and 
can easily outshine it. This holds true in all Sommerfeld models that we have considered. 
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Figure S6: The effect of substructure on the local annihilation rate. A: The mean and maximum number of Via 
Lactea II subhalos inside 100 randomly placed spheres 8 kpc from the halo center versus the radius of these sample 
spheres. The mean subhalo occupancy becomes unity at r = 3 kpc. B: The density "clumping factor" (p 2 )/(p) 2 
over the 100 sample spheres. C: The ratio of the subhalo to host halo contributions to the annihilation luminosity 
for three representative Sommerfeld models (none, 1/v, and 1/v 2 ). Only subhalos resolved in our simulation are 
accounted for. 

S6 Local Luminosity Boost From Substructure 

Here we assess the role that nearby subhalos play for the Sommerfeld-enhanced production 
of high energy electrons and positrons. Because these energetic particles lose energy as they 
diffuse through the Galactic magnetic field, only those that are produced within a few kpc of 
Earth are of interest. 

We considered 100 spheres of radius 20 kpc, with randomly positioned centers 8 kpc from 
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the host halo center. For each of these spheres we determined the cumulative number of subha- 
los N suh , a "clumping factor", defined as (p 2 ) / (p) 2 , and the ratio of total subhalo to host halo 
luminosity as a function of enclosed radius in the sphere. For the subhalo luminosity we used 
the analytical NFW estimate, as explained in Section [SH Figure [S6] shows the mean and maxi- 
mum values of these quantities over all 100 sample spheres. The low local subhalo abundance 
is reflected in a small mean iV sub within a few kpc of the Sun. Only 3 of the 100 sample spheres 
have any subhalos within 1 kpc of their center. The mean subhalo occupancy becomes unity at 
3 kpc, but one has to go out to 7 kpc before every single sphere contains at least one subhalo. 
The clumping factor captures the enhancement of the annihilation luminosity compared to a 
homogeneous density background. It has contributions from the overall density stratification, 
from the Poisson noise of the density estimator, and from unbound and bound substructure. 
The sharp rise towards 8kpc is due to the cuspy nature of the host halo density profile. The 
bottom panel of Figure [S6] shows that without Sommerfeld enhancements the subhalos do not 
contribute significantly to the local annihilation luminosity. In a typical S ~ 1/v model, how- 
ever, subhalos contribute on average about half as much as the host halo, and in rare cases 5 
times more. For models on resonance (S ~ 1/f 2 ), the subhalos completely dominate the host 
halo, and provide on average 20 times as much luminosity as the host halo. Again we have 
neglected the possible additional contribution from subhalos below our simulation's resolution 
limit. 
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